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ABSTRACT 

We present an idealized, semi-empirical model for the evolution of gravita- 
tionally contracting molecular clouds (MCs) and their star formation rate (SFR) 
and efficiency (SFE). The model assumes that the instantaneous SFR is given 
by the mass above a certain density threshold divided by its free-fall time. The 
instantaneous number of massive stars is computed assuming a Kroupa IMF. 
These stars feed back on the cloud through ionizing radiation, eroding it. The 
main controlling parameter of the evolution turns out to be the maximum cloud 
mass, M max . This allows us to compare various properties of the model clouds 
against their observational counterparts. A giant molecular cloud (GMC) model 
(M max ~ 10 5 M Q ) adheres very well to the evolutionary scenario recently in- 
ferred by Kawamura et al. (2009) for GMCs in the Large Magellanic Cloud. A 
model cloud with M max 2000 M Q evolves in the Kennicutt-Schmidt diagram 
first passing through the locus of typical low- to-intermediate mass star-forming 
clouds, and then moving towards the locus of high-mass star-forming ones over 
the course of ~ 10 Myr. Also, the stellar age histograms for this cloud a few Myr 
before its destruction agree very well with those observed in the p-Oph stellar 
association, whose parent cloud has a similar mass, and imply that the SFR of 
the clouds increases with time. Our model thus agrees well with various observed 
properties of star-forming MCs, suggesting that the scenario of gravitationally 
collapsing MCs, with their SFR regulated by stellar feedback, is entirely feasible 
and in agreement with key observed properties of molecular clouds. 

Subject headings: ISM: clouds — ISM: evolution — Stars: formation 
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Introduction 



A crucial ingredient in understanding the star formation efficiency (SFE) of giant 
molecular clouds (GMC) is the study of their evolution, from their formation to their 
destruction by the massive stars they form. A still unsolved problem is whether the GMCs 
are in approximate virial equilibrium, or rather they are in gravitational contraction. 
Initially, Goldreich & Kwan (1974) proposed that the supersonic linewidths observed in 
GMCs correspond to global gravitational contraction, but Zuckerman & Palmer (1974) 
readily argued that if all the molecular gas in the Galaxy were in free fall, then the 
total star formation rate (SFR) in the Galaxy would be about two orders of magnitude 
higher than observed (we will refer to this as the "SFR conundrum"). Zuckerman & 
Evans (1974) subsequently suggested that the linewidths could correspond to small-scale 
turbulent motions, giving rise to the notion that clouds are q uasi-equilibrium entities, a 
notion that has survived until today (see e.g., the reviews by 



Mac Low fc Klessen 



McKee Ostriker 



2007J). 



2004; 



Since then, most theoretical models of star formation (SF) have been based on the 
assumption that turbulence provides support against the clouds' self-gravity, and allows 



maintaining a low g 


lobal SFR (e.g., 


Norman & Silk 


1980 




McKee 


1989 




Matzner 


ZOoi 


Krumholz & McKee 


2005; 


Li & Nakamura 


2006; 


Nakamura & Li 


2007 




Wang et al 


2010) 



In these models, the turbulence is maintained by the energy feedback into the cloud from 
the star s it forms. One i nteres ting model where strict equilibrium was not assumed, was 
that by 



Krumholz et al. 



(120061 ) . where the fully time-dependent Virial Theorem was solved 
numerically for a cloud under the influence of its self-gravity and the pressure produced by 
feedback from H II regions, although the cloud was restricted to have a spherical geometry, 



and mass loss by the cloud due to ionizing radiation by massive stars was not considered. 
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Those authors found that clouds undergo a few expansion-contraction oscillations, until 
they are finally dispersed, and the SFEs over the clouds' lifetimes were found to be ~ 5 
- 10%. More recently, a similar model, with the s ame restrictions but in cluding mass 
accretion from the environment, was considered by 



Goldbaum et al. 



(1201 1( 1 . In this model, 



the clouds again reach virial equilibrium, and maintain roughly constant column densities. 

However, recent theoretical and observational evidence has suggested a re- 
turn to the global gravitational contraction scenario of Goldreich &: Kwan (1974 ) 



Har tmann. Ballest eros-Paredes. Bergin 



2001 



Burkert fc Hartmann 



2004 



Hartmann Burkert 



2009 



2007 



Peretto, 



Vazquez- Semadeni et al. 



2007, 



iennebelle Andre 



2009 



Csengeri et al. 



2007 



Galvan- Madrid et al 



2010 



Schneider et al. 



2010|). 



Additionally, IClark fc Bonnelll (j2005[ ) have suggested that molecular cloud turbulence, 
rather than directly producing Jeans-unstable clumps, only produces the seed nonlinear 
density fluctuations for subsequent gravitational fragmentation, which proceeds on different 
timescales due to the spatia l variations on the local free-fall time induced by the turbulence 



( iHeitsch fc Hartmann 



20081 ). But then, if we again allow for clouds and their substructures 
to be in gravitational contraction, it is necessary to find a solution for the SFR conundrum 
in this scenario. The semi-empirical model presented here investigates whether stellar 
feedback can accomplish this. 



Our model is motivated by the numerical simulations of 



Vazquez-Semadeni et al 



( 120101 ) . who have investigated the evolution of clouds formed by the collision of warm neutral 
medium (WNM) cylindrical streams, including stellar feedback from ionization heating 
from massive stars. Those authors found that the clouds are in general not stabilized by the 
feedback, but rather are either dispersed or continue to contract globally, depending on their 
mass. Here, we construct a model that attempts to capture the phenomenology observed 
in those simulations, scanning the space of the parameters that determine the physical 
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properties of the cloud. iRosas-Guevara et al.l (120101 ) have also prese nted a parameter -space 



Dibetal 



(1201 If ) have 



study, but using numerical simulations, and without feedback, while 
produced an analytical model similar to ours, but aimed at investigating the effect of 
varying metallicity on the SFE. 

)e the general model, wh ich 



Vazquez-Semadeni et al. 



The plan of the paper is as follows. In section §2] we descri 
we th en calibrate against a fiducial numerical simulation from 
( 120101 ) in £j3j In £j5]we compare the calibrated model against various observational properties 
of both large and small MCs, parameterized only by their mass. In §|6] we present a 
discussion, and finally a summary and our conclusions in §0 



2. The Model 



We construct a model for studying the SFE in a thin cylindrical cloud undergoing 



gravi t ational contraction, as observed in various simulations ( TVazquez-Semadeni et al 



2010 



2011 



20071 . 



Heitsch fc Hartmannl 120081 ). The system is schematically illustrated in Fig. 



[TJ In the simulations, the collision of warm neutral medium strea ms nonlinearly triggers 



thermal instability, forming a thin cloud of cold atomic gas (e.g. 



1999 



Koyama &: Inutsuka 



Hennebelle fc Perault 



2000 



2002 



Walder fc Folini 



the combined action of various dynamical inst abilities ( 



Koyama &: Inutsuka 



2002 



Heitsch et al 



2005 



2000 f), which 



Decomes turbu 



Hunter et al. 



1986 



Vazquez-Semadeni et al. 



Vishniac 



ent by 



1994 : 



20061 ). The cloud 



soon begins to contract gravitationally as a whole. However, before this global collapse is 
completed, some local, nonlinear (i.e., large-amplitude) density enhancements produced 
by the initial turbulence manage to collapse on thei r own, since their local free-fall time 



is sho rter than the average one for the entire cloud (IHeitsch fc Hartmann 



2008 



Pon et al. 



20111 ). These local collapses thus involve only a fraction of the cloud's total mass. Also, 



we assume that the newly-formed stars feed energy back into the cloud. We only consider 
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the ionizing radiation from massive stars, since this radiati on is probably the dominant 



mechanism of stellar energy injection at the scale of GMCs (IMatzner 



20021 ) 




Fig. 1. — Model setup. Cylindrical streams of WNM are assumed to collide head-on to 
form first a flattened CNM cloud, which proceeds to collapse, becoming molecular and star- 
forming in the process. The main parameters of the model are the radius of the cylinders 
Ri n f , and the density and velocity of the inflowing warm gas. 



In what follows, we investigate the competition between the cloud's gravitational 
contraction and its destruction by the mass consumption by star formation (SF) as well as 
by the ionization produced by the newly formed massive stars. Below we describe how we 
calculate the contributions from these processes. It is important to note that we do not 
follow the chemistry, but rather consider that all the cold gas, either atomic or molecular, 
is involved in the gravitational contraction and, eventually, star formation. 



2.1. Mass accretion 

In our scenario, the cloud's mass (Mq) at time t is given by 

M c (t)= [ M M (t') dt'-Ms(t)-Mi(t), (1) 
J o 
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where M- m f(t) is the mass accretion rate onto the cloud from the WNM inflows, Ms(t) is the 
total mass in stars, and M\(t) is the mass ionized by stellar feedback. Note that all these 
masses are considered to be functions of time. We assume that 

M inf (t) =2 Phli v ml (7rR 2 c ), (2) 

where v in f is the inflow speed, p in f = n in f /in m H is the inflow mass density (with n inf the 
number density of the inflows, /in the mean atomic weight of the diffuse gas, and hih 
the atomic hydrogen mass). Note that it Rq is also the cross-sectional area of the thin, 
cylindrical, cold and dense cloud that forms by thermal instability at the layer compressed 
by the inflows. We take the cloud's radius Rc(t) as being initially equal to the inflow radius 
R m {, and to later decrease as the cloud contracts gravitationally. The factor of 2 in eq. (j2J) 
represents the fact that there are two inflows, one on each side of the forming cloud. 



2.2. Mass in Stars 

We assume that the SFR is given by the ratio of the gas mass in the high-density tail 
in > n SF ) of the density distribution produced by the turbulence in the cloud, to its local 
free-fall time, tff(nsF) = a/37t /32G/imH^sF- We refer to hsf as the threshold density for 
star formation, and denote by /(f) the fraction of the cloud's mass that is at densities above 
n SF (discussed in §2.6p . We then have 

SFR(t) = i^*) /(f), (3) 
Wsf) 

so that the mass in stars at time f is 

M s (f) = jf SFR(f')df = jf* f(t') df. (4) 

We assume that the initial density of the cloud is that of the cold neutral medium 
(CNM), in balance with the sum of the thermal and ram pressures of the inflows, as 
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described in 



Vazquez-Semadeni et al 



( 120061 ). Typically, t^cnm ~ 100 cm 3 , and we use 



a mean molecular weight of /i = 2.35, adequate for molecular gas. As the cloud evolves 
by contraction, the mean density evolves as determined by its mass and size, so that 
p = Mc/nRch, where h is the cloud thickness (cf. §2.4[) . 



2.3. Ionized Mass 



To model the cloud evaporation by massive stars, we use the results from 



Franco. Shore fc Tenorio-Tagld ( 119941 ) . These authors found that the cloud evaporation 



rate by a massive star near the cloud surface is 



M I)Sur (t) « 2tcRI m H nc Si i ( 1 + 



5c s ,ity/5 

2-Rs,o 



(5) 



where t is the age of the massive star, c Sj i is the sound speed in the ionized gas, n — pj pmu 
is the mean number density of the molecular cloud (MC), and i?s,o is the initial Stromgren 
radius of the massive star in the cloud (Rs,o = \ 3S* / 4:%ap.(2n) 2 } 1 / 3 , with S± a represen tative 



value of the UV Lyman-continuum photon flux ([Franco. Shore fe Tenorio-Tagle 



199J), and 



ctB the recombination coefficient for the ionized gas). We assume that i?s,o is reached 
immediately at the instantaneous mean density of the cloud, which, in our model, is 
continually increasing as the cloud contracts. Therefore, over a short time interval At 
(between t and t + At), over which the cloud's density can be assumed constant, a massive 
star near the cloud's surface can ionize a mass AMi sur = Mj jSW (t)At. 

In addition, we consider that, during the time interval At, the cloud forms 
A-/Vob(£) = £obSFR(£) At/ (Mqb) new massive stars of average mass (Mob), wh ere xqr 



i s the mass fraction of massive stars, which we calculate assuming an IMF from 



Kroupa 



(120011 ). with lower and upper mass limits of 0.01 and 60 M respectively. Defining a star as 
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"massive" if it has a mass M > 8M , we obtain xob = 0.12 and a mean massive-star mass 
(M OB ) = 17 M . Furthermore, in eq. © we take S, = 2 x 10 48 s-\ c s>1 = 12.8 kms" 1 , and 
qb = 2.6 x 10 -13 cm -3 s -1 where we have assumed that the temperature of ionized gas in 
the HII region is 10 4 K. 

With the above considerations, and discretizing the time variable, the ionized mass at 
time ti is given by 

i 

Mj(U) = Mife-i) + A#oB(*i)AMi, Bur (* i ), (6) 

j=k 

where is the time at which the oldest remaining OB stars were formed, k = i— int (£ob/A£), 
'int' is the integer function, and toB = 5 Myr is the main-sequence lifetime of our 
representative OB star. The second term in the right-hand side of eq. (j6]) thus gives the 
mass ionized over the time interval At = ti — by the stars formed between and the 
present time, ti. Note that we have neglected the possibility that part of the ionized gas 
can recombine and return to the cloud, as we expect the ionized gas to escape to the diffuse 
medium. 

Finally, to account for the fact that the massive stars are born in dense environments, 
in eq. fl5]) we take the density as the maximum between the instantaneous cloud density 
n(t) and 10 3 cm -3 (typical of clumps), in order to avoid an over-ionization when the cloud 
density is low. 



2.4. Global gravitational collapse. 



Following the trend seen in the numerical simulations, we assume that our clouds evolve 
in two stages. First, a mass-growth stage occurs, during which the cloud (initially of zero 



mass) increases its mass by accretion 
only its thickness increases (see, e.g., 



rom the WNM at constant radius and density, so tha ; 



Vazquez-Semadeni et al. 



2006 



2007 



Folini fc Walder 
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20061 ). until it reaches its thermal Jeans mass. At that point, the second stage begins 
during which the cloud undergoes global gravitational contraction. For circular modes in 
a sel: 



1985|): 



-gravitating isothermal sheet of finite thickness, the Jeans mass is given by (1L arson 



Mj = 4.67 



G 2 S 



(7) 



where c s is the cloud sound speed, and £ = M(t)/nR 2 (t) is the surface density, with M{t) 
and R(t) being the instantaneous mass and radius of the cloud, respectively. In order to 
calculate the sound speed we first compute the cloud's temperature. To do this, we use the 
fit by Koyama & Inutsuka (2002; see also th e note in Vazquez- Semaden i et al. 2007) to the 



heating and cooling processes considered by 



Koyama &: Inutsuka! (120001 ) . This allows us to 



solve for the temperature of thermal equilibrium (heating = cooling) as a function of the 
density. In this way, we get temperatures of T rj 42 K for a density of n = 10 2 cm -3 and 
T ps 7 K for n ~ 10 7 cm" 3 . 

During the mass-growth stage, over which the cloud's radius remains constant 
(R = Rinf), the cloud's thickness h is given by the condition of constant number density: 

Once the collapse begins, we assume that the thickness remains constant at the final 
value achieved during the growth stage, and that its volume density increases only 
due to the radial contraction, as suggested by simulations including self-gravity (e.g., 



Vazquez-Semadeni et al. 



2007 



2010 



2011 



Heitsch fc Hartmann 



20081 ). This assumption 



is equivalent to assuming that the average thickness of the cloud is much smaller than its 
Jeans length throughout its evolution. In general, this is a good approximation until when 
the cloud has contracted to radii of a few pc. 

To determine the radial evolution during the contraction stage, we first calculate the 
acceleration at the cloud's edge. We take a reference frame with its origin at the cloud's 
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center and with its a;-axis along the inflow direction (i.e., perpendicular to the plane of 
our flattened cloud). We integrate over mass elements p dzdydx located at a distance 
[(R — z) 2 + y 2 + x 2 } 1 / 2 from the edge. At a certain time £j, the acceleration at the cloud 
edge is 

rRi ry/R?-z* rh/2 R _ z 

a(U) = 2 G ft * I A, m _^ + yl + xY/1 dx (9) 

where pi = p(ti) and R4 = Rq^). We solve the first integral analytically, while the second 
and third ones are solved numerically by the composite Simpson rule. After a small time 
increment At = U — the change in the cloud radius is: 

R l+1 =Ri- v ^ At - ^a{Ri)At 2 (10) 

where the instantaneous velocity at time tj is v 0>i = J2] =0 ci(Rj)At for constant At. 

Finally, we introduce a correction factor, representative of the fact that true 
gravitational collapse of a gaseous mass does not occur in strict free-fall, since t hermal 



pressu re is never completely negligible, as pointed out in the pioneering work by 



Larson 



(11969 



Appendix C). There, he reported that the collapse of his simulations occurred in a 
time longer than the free-fall time by a factor of 1.58. Thus, at each timestep, we divide the 
radius given by eq. ( TTOj) by a factor /l, which we refer to as the "Larson parameter", and 
calibrate in 33l 



2.5. Density distribution 

A key ingredient in our evolutionary model is the fraction of dense gas that is 
participating in the SF process, and therefore so is the evolution of the probability density 
function (PDF) of the cloud's density field. 

In our model, we consider that our clouds are born transonically turbulent, as a 
consequence of the various instabilities at play in the compressed layer between the streams 



(IVishniad 11994 : IWalder fc Folinil 120001 : IHeitsch et al.N2005L 120061 : IVazquez-Semadeni et al 
20061 ) . Also, in the absence of direct stellar irradiation, the temperature in the cold gas 
varies at most by factors of a few for densities 100 < n < 10 7 cm -3 , and thus, as a first 
approximation, we consider it to be isothermal. Therefore, we assume that the density 
field within the cloud is initially charact erized by a lognormal PDF, appr opriate for 



supersonically turbulent, isothermal gas (IPassot fc Vazquez-Semadeni 
then given by 



Pis) 



1 



exp 



S - S T 



2al 



where 



ln(p/p), 



HPp/p) = "-'A 



19981 ) . The PDF is 
(11) 
(12) 



with p p the peak density, and 



a 2 s = ln(l + b 2 M 



2 ) 

rms/i 



(13) 



where h is a proportionality constant related with the compressibility induced by the 



turbu 



1994 



2008|). 



ent forcing, which 



Padoan et al. 



1997 



'or simplicity we take equal to uni t y (see e.g., 



Passot fc Vazquez-Semadeni 



1998 



Vazquez-Semadeni 



Federrath. Klessen &: Schmidt 



However, more recent numerical and observational studies suggest that the density 



the collapse, but rat 


ler develops a power-law tail at high densities 


Klessen 


200(J: 


Dib & Burkert 


2005: 


Vazauez-Semadeni et al. 


2008: Kainulainen et a . 


20oJ: 


Kritsuk et al. 


2011; 


Ballesteros-Paredes et al. 


2011b 


). In particular, 


Kritsuk et al. 


(2 


011 


) have suggested 



that the final slope should be in the range [3/2,7/4], but at the present time we know 
of no theoretical prediction as to how the slope nor the transition point between the 
lognormal and the power law should evolve in time. We experimented with various options 
for modeling the evolution of the PDF's power-law tail, but found the behavior to be 



13 



very sensitive to the parame t ers us ed, while we had no physical ways of constraining 



them. Finally, 



Kritsuk ct al. 



(1201 if ) proposed that the origin of the power-l aw tail is the 



development o f lo cal collaps ing flows, that may have either Larson-Penston (ILarson 



Penston 



1969 



19691 ) or IShul (119771 ) density profiles. In this case, the power-law tail is an effect 
of the gravitational collapse rather than its cause, and thus it should not be counted as 
providing turbulent seeds for future collapses. For all of these reasons, we do not consider 
the power-law form of the PDF, and stick to the lognormal. 

It is worth noting that the density PDF we consider here refers only to the cold 
(approximately isothermal) gas that makes up the cloud, and not to the entire gas contents 
of the system. This means that this PDF is not directly comparable to that observed in 
the numerical simulations of the same process, which corresponds to thermally-bistable 
gas. Thus, we did not consider the possibility of using a PDF extracted from numerical 
simulations, either. 

We model the evolution of the lognormal density PDF as follows. First, as indicated 
by eq. ( TT2"j) . the mean of the PDF is given by the cloud's mean density, given by 
p = Mc/tiRqIi. In turn, we prescribe that p varies as follows. During the mass-growth 
stage, it remains constant, at 100 cm -3 . Once the contraction stage begins, it increases, 
causing the density PDF to shift to higher values. Eventually, however, the SFR becomes 
large enough that the ionization by newly born massive stars reduces the cloud's mass 
rapidly enough as to cause the PDF to shift back towards lower densities again. 

Second, the standard deviation of the PDF is determined by the turbulent rms Mach 
number, as indicated by eq. ( 1131) . Unfortunately, the evolution of the turbulent component 
of the Mach number remains rather uncertain. Standard relations, such as Larson's 
(1981) velocity dispersion size cannot be assumed here. Indeed, in our collapsing-cloud 
scenario, the majority of the velocity dispersion is due to the contracting motions 
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( IBallesteros-Paredes et al.ll201ll ) rather than to random turbulent motions, and should 
not be counted as turbulence capable of producing new density fluctuations suceptible of 
subsequent collapse. 

Of course, it is natural to assume that a fraction of the kinetic energy in the collapsing 



motions will be transferred to random motions, but this prob 



in the case of gaseous media. 



Vazquez-Semadeni et al 



em remains largely unexplored 



( 119981 ) numerically investigated the 



scaling of the non-collapsing component of the velocity dispersion in collapsing spherical 
clouds. Those authors found that the turbulent velocity dispersion scaled as p x , with 
x G [1/4, 1/2], depending on the particular setup of the collapse and on the presence of 
magnetic fields. However, that study was performed at low resolutio n and was restricted to 



spheri cal geometry, so it cannot be taken as definitive. More recently, 



Klessen fc Hennebelle 



(120 101 ) suggested that the rate of kinetic energy injection by the warm neutral streams 
feeding a molecular cloud is at least one order of magnitude larger than the rate of 
turbulent energy dissipation within the clouds. However, these estimates are not enough 
to properly constrain the evolution of the turbulent kinetic energy in our model clouds. 
Thus, we simply take a constant value of the turbulent Mach number, which represents 
a compromise between turbulent decay by dissipation and feeding of the turbulence by 
transfer from the collapsing motions. This thus implies that the width of the PDF in 
our model remains constant through the cloud's evolution as well. We take this constant 



value of t 



re turbulent Mach number as Ai 



medium (IHeiles fc Troland 



3, the canonical value for the cold neutral 



2OQ3J). 



2.6. Star-forming mass fraction 

Finally, as mentioned in §2.2[ we assume that only gas with number density higher 
than a threshold value tt-sf participates in the instantaneous SF process. This amounts to 
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assuming that usf is sufficiently larger than the cloud's mean density n as to guarantee 
that tff(rigF) <C tff(n), where tg(n) is the free-fall time corresponding to density n. In a 



sense, this prescriptio n may be t 



numerical simulation ( Bate et al. 



lought of as the mode 



1995 



Federrath et al. 



s an alogue of sink particles in a 



20101 ) . in which gas at sufficiently 



high densities is replaced by point mass particles representing collapsed objects. 

Thus, for the assumed lognormal PDF the fraction of dense gas that forms stars is 
given by 



/(*) 



erf 



2s SF (t) -a\ 



V2a s 



(14) 



where srf 



Elme green 



mfisp/ff), and pm? is the volume density 



2002 



Krumholz fc McKee 



2005 



Dib et al 



corre sponding to n SF (see also 



201 111. The value of wqf is 



calibrated against a numerical simulation by IVazquez-Semadeni et al 



d2010h in 



2.7. Temporal evolution 

With all the above ingredients, we discretize eq. ([I]) as 
MdU) = Y,M in{ (t 3 ) A t-J2 YW) [ M ^-^ + J2 ^osfeOAMi^fa)] (15) 

3=0 j=0 S{ - j > j=k 

We integrate this equation numerically over time, taking At = 4 x 10~ 4 Myr. We choose 
this value as a reasonable compromise between speed and accuracy, after experimenting 
with timesteps down to 0.01 times this value, from which we found (with arbitrary free 
parameters) that the chosen value gives already well-converged values of the final cloud 
masses, the variable that turned out to be most sensitive to the size of the timesteps. 

We consider that the cloud's evolution ends either because the cloud gets completely 
evaporated by the massive-star ionization or because its average density reaches n^p, a 
point at which the remaining mass is fully converted into stars and the cloud disappears. 
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Calibration of the model 



With all the model ingredients de fined, we now proceed to calibra te it by matching it 



to one of the numerical simulations by 



Vazquez-Semadeni et al 



(12010h . We consider the 



simulation labeled SAF1 in that paper, which most resembles the physical conditions built 
into our model. The label SAF indicated that it contained small- amplitude initial velocity 
fluctuations (~ 2% of tw), which were sufficient to trigger the instabilities responsible for 
turbulence production in the forming cloud, but not large enough to significantly distort its 
sheet-like geometry. 

The simulation considered the collision of two cylindrical WNM streams with 
n in f = lcm~ 3 , TLf = 5000 K, inflow radius of i?j n f = 64 pc, and an initial inflow speed 
v- m { = 7.5 km s 1 J^J Thus, for our calibration purposes, we use these same values for our 
model, except for t> in f. We do this because, in the simulation, the inflow speed decreases 
over time, because the rear end of the inflowing cylinder leaves a vacuum behind it. As a 
consequence, in the reference frame of the inflow, the inflow's rear end tends to re-expand 
into this vacuum, producing a velocity gradient along the cylinder. Instead, in the model 
we use a constant inflow speed, and thus we take a smaller value, representative of the 
time-averaged speed in the simulation. We find that a value of v- m f = 4.5 km s -1 provides 
the best match for the evolution of the simulation cloud's mass. Also, we find that 
the temporal evolution of the cloud's radius is best matched with a value of the L arson 



param eter of /l = 1.7, very close to the original value of 1.58 originally proposed by 



Larson 



( 119691 ). We show the evolution of these quantities for both the simulation and the model in 



Fig. |2j Note that the evolution for the model with the ionizing stellar feedback turned both 
on and off is shown. We refer to the model with feedback off as the LN-FO case, and to the 



1 Note that there are typographical errors in the numbers reported by 



Vazquez-Semadeni et al 



( 120 lOh . The values given here are the correct ones. 
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one with the feedback on, as the LN-F1 case. 
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Fig. 2. — Left: Evolution of the dense gas mass, ionized gas mass, and mass in stars for the 
model cloud with parameters that best match the corresponding quantities (except for the 
i onize d mass, which is not measured) in the SAF1 simulation of IVazquez-Semadeni et al. 
( I2010l ). The green lines represent the simulation, the blue lines represent the model cloud 
with stellar feedback turned off, and the black lines represent the model cloud with the 
stellar feedback turned on, with usf = 10 6 cm" 3 . The vertical line represents the beginning 
of global collapse. Right: Evolution of the cloud radius in the model and in the simulation, 
with the same color coding as in the left panel. 



Next, keeping the above parameters fixed, we vary the n^F parameter, the density 
threshold for star formation, to match the SFR and SFE of the simulation, with the SFE 
being defined as 

SFE(t) = -— - — - , (16) 



M c (t) + M s (t)+Mi(ty 
where Mq is the dense gas mass, Mg is the total mass in stars, Mi is the total mass ionized 

by stars, all quantities being time-dependent. Figure |3] shows the evolution of the SFR and 

SFE of the model for ngp = 10 5 , 10 6 , and 10 7 cm -3 , and compares them with the evolution 

of the corresponding quantities in the simulation, showing that the best match is obtained 
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with a value risF = 10 6 cm 3 , which we use in the rest of the paper. This value is reassuring 
since, on the one hand, it i s comparable to the sink-formation d ensity threshold used in 



the numerical simulations (jVazquez-Semadeni et al. 



2007 



20101 ). and on the other, it is 



high enough that the material above those densit i es can be safely assumed t o be locally 



gravitationally bound ( iGalvan-Madrid et al. 



2007 



Heitsch fc Hartmann 



2008|). 



LN — F1 model with n 

o 4 |- i 5 cm" 

s cm"' 
1 7 cm"' 




Time [Myr] 



Time [Myr] 



Fig. 3. — Left: Evolution of the star formation rate (SFR) for the model cloud for the cases 
rig g = 10 5 , 10 6 , and 10 7 cm" 3 , com pared to the evolution of the SFR in the SAF1 simulation 



of 



Vazquez-Semadeni et al. 



( 120101 ). Right: Same as the left panel, but for the star formation 



efficiency (SFE). The green lines represent the numerical simulation. 



4. Discussion of the calibrated model's evolution 

Once the model parameters have been calibrated with the SAF1 simulation, it is 
illustrative to discuss the general features of the model cloud's evolution, which hold 
qualitatively for the other cases we explore in £JS1 where varying only the inflow radius 
we are able to match and explain several observed features of molecular clouds and their 
star-forming activity. 
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First, we note from the right panel of Fig. [2jthat the radius of the model with feedback 
(LN-F1) evolves slightly more slowly than that of the case without feedback (LN-FO). This 
is because the stellar feedback erodes the cloud through ionization, thus reducing its mass, 
which in turns causes a lower gravitational acceleration, thus slowing the collapse. We also 
note that the radius of the LN-FO case approaches zero at late times, while that of the 
LN-F1 ends at a finite radius, implying that the evolution is terminated because the cloud 
is completely evaporated before it reaches zero radius (see below). 

Second, from the left panel of Fig. [5J we note that the dense gas mass of both the 
LN-FO and LN-F1 models decreases at late times. This is because even in the non-feedback 
LN-FO model, gas is consumed due to the conversion of gas to stars. However, we see that 
the mass consumption is much more abrupt in the LN-FO model, while in the LN-F1 model 
it proceeds more slowly. This indicates that the feedback inhibits star formation to the level 
that the mass consumed by ionization from the feedback is less than the mass that would 
be consumed by star formation were there no feedback. Nevertheless, note that the dense 
gas mass in the LN-F1 model approaches zero at the end of the evolution, indicating that 
all of the cloud's mass is used up by the combined action of ionization and star formation. 



This corresponds to the well known fact that clusters eventually 



destroy their natal c 



19891). 



oud, 



and are left with no gas around them after several million years (iLeisawitz et al.l 

Finally, note from Fig. H] that both the SFR and the SFE increase at an ever faster 
pace until the end of the simulation in the non-feedback case LN-FO, while their growth 
slows down in the case of LN-F1. In this case, the SFR eventually begins to decrease and 
goes to zero at the end of the evolution, as the dense gas mass is completely consumed by 
the ionization, leaving no further fuel for star formation. Instead, the SFR in the LN-FO 
case would reach a singularity at a finite time were it not for the time discretization of our 
model. 
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Fig. 4. — Evolution of the SFR (left) and the SFE (right) of the calibrated model with 
n SF = 10 6 cm" 3 for with feedback turned off (LN-FO, blue line) and cL CcLSG with 

feedback on (LN-F1, black line). The green lines represent the numerical simulation. The 
LN-FO model has accelerating SFR and SFE, while they decelerate for model LN-F1. See 
also discussion in §6.3! 

5. Comparison with observations and previous work 

We now proceed to compare the results of our model with related observational 
results. For the various comparisons, we vary only the inflow radius parameter, -Rj n f, which 
determines the maximum dense gas mass attained by the model. 



5.1. Evolutionary stages 



As a first ca s e in p oint, we consider the evolutionary stages of the clouds. Recently, 



Kawamura et al. 



( 120091 ) have suggested that GMCs in the Large Magellanic Cloud undergo 



four evolutionary stages. In the first stage (Type I, with a duration of ~ 7 Myr and a 
median mass M t i ~ 10 4 ' 8 M ), the GMCs show no sign of massive star formation. In the 
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second stage (Type II; 14 Myr; M tII ~ 10 5 2 M ), the GMCs have only HII regions, while 
in the third (Type III; 6 Myr; M tm ~ 10 5 - 4 M ), the GMCs contain both HII regions and 
clusters. Finally, the last stage (IV) corresponds to the time when the GMCs have been 
completely dispersed, and only young clusters and/or SNRs are found. 

Because the masses of the clouds in their sample are near 10 5 M Q , we choose an inflow 
radius i?; n f = 100 pc, which gives a maximum cloud mass (M max ) of slightly over 10 5 M , 
thus making it directly comparable to their cloud sample. Figure |5] shows the evolution 
of both the number of massive stars in the model cloud (top panel) and the masses of the 



dense gas, ionized gas, and stell ar components o 



with the cloud types defined by 



Kawamura et al. 



the c 



oud (bottom panel). For comparison 



(120091 ). here we define Type I as the epoch 



when the model cloud has less than one massive star, Type II as the period when the 
cloud has less than 20 massive stars, and Type III as the period when the cloud has more 
than 20 massive stars. We see that the model cloud spends ~ 5 Myr as a Type I, ~ 12 
Myr as a Type II, and ~ 10 Myr as a Type III, noting that after such a time the cloud's 
mass has decreased by more than a factor of 2, and may be considered to be on its way to 
disappearing. 

We also note from the bottom panel of Fig. [5j that the dense gas mass of the model 
cloud varies only moderately during the time it spends as an either Type I, II or III cloud, 
although a net increase from Type I to Type II is apparent. Moreover, because the cloud's 
mass decreases by over a factor of 2 while it is in the Type III stage, a si gnificant scatter in 
cloud masses is expected in this class, as is indeed observed in Fig. 12 of 



Kawamura et al 



fl2009h . 



From the above discussion, we thus conclude that the evolution of our model GMC, 
with i? in f = 100 pc, compares well, both qualitat i vely a nd quantitatively, with the 



evolutionary scheme proposed by 



Kawamura et al. 



(120091 ) for GMCs in the LMC. 
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Fig. 5. — Top: Evolution of the number of massive stars in the model cloud with 
R\- nf = 100 pc, showing th e periods which roughly correspond to the cloud types defined 



by 



Kawamura et al 



( 120091 ). Bottom: Evolution of the dense gas, ionized and stellar masses 



for this cloud. 
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5.2. Fiducial model vs. OMC-1 



Next, we compare our model cloud with the physical conditions of real star-forming 
regions. As an example, we choose the clump known as OMC-1 in the Orion Molecular 
Cloud, which can be considered a typical massive star-forming re gion, and has be en 



extensively studied. The total gas mass in OMC-1 is ~ 2200 M (jBally et al.l 119871 1. the 
size is ~ 1.35 pc (which implies a number density ~ 1.54 x 10 4 cm~ 3 ), and the mass in stars 
is ~ 500 M , which implies that the average SFR is (SFR) ^ 2.5 x lO _4 M yr _1 , assuming 



a stellar age spread of < 2 Myr (see 



Vazquez- Semadeni et al. 



2009 



and references therein). 



- <SFR> over At = t-2Myr 
LN-FO mode 
LN-F 




R [pc] 

Fig. 6. — Evolution of a model cloud with i? in f = 10 pc in a SFR vs. size (R) diagram. 
Note that the size axis is reversed, so that the cloud evolves from left to right as it contracts 
gravitationally. Also plotted is the locus of the OMC-1 star-forming region in this diagram, 
showing that the model comes very close to that locus towards the end of its evolution. 

To compare with this, we choose a value of i?i n f = 10 pc for our model, which at 
t = 27.6 Myr has the same density as OMC-1, making it directly comparable to the latter. 



At this time, the model cloud has a mass M ~ 1800 M , a mass in stars ~ 200 Mm, and size 



L 0, 
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(diameter) « 1.9 pc. Moreover, the average SFR over the last 2 Myr (the age dispersion 
used to compute the SFR of OMC-1) is (SFR) w 100 M Myr^ 1 (see Figure ED- These 
values are in very good agreement with the estimates for OMC-1, suggesting that our 
evolutionary model correctly describes this cloud. 



5.3. Kennicutt-Schmidt relation. 



Another possible point of comparison of our model with observational dat a is provid ed 



by the so-called Kennicutt-Schmidt relation. Ever since the seminal paper by 



Schmidt 



(119591 ). it has been well known that there exists a relationship between the SFR and the 



gas density in galaxies. Four d ecades late r, collec ting data from various surveys of nearby 



normal and starburst galaxies, 



Kennicutt 



(ll998bl ) found a clear correlation between the 



galaxy-averaged SFR surface density (Ssfr) and the galaxy-averaged total gas surface 
density (S gas = £ H i + S H2 , where S H i and S H2 are the HI and H 2 surface densities, 
respectively) of the form S SFR oc £^ s , with iV rs 1.4. 

Recent observations of external galaxies with high spatial resolution (on scales of 
fractions of kpc) show th at Srfr scales a l most linearly with Sh 2 , while no clear correlation 



exists with Shi (see e.g., 



individual clouds (e.g., 



Wu et al. 



Evans et al. 



2005 



2009; 



Big iel et al. 



Heiderman et al 



values of Srfr t han those implied by the fits by 



Kennicutt 



Wu et al. 



(120051 ). Moreover, the SFRs derived by 



2008). However, observations of 



2010) sy stematically show larger 



1998b), 



Heiderman et al. 



Big iel et al. 

"IT 



(J2008|), and 



2010) for their massive 



clump sample were obtained using extragalactic methods (taken from IWu et al.l 120101 ) 



and they warn that this could cause the SFRs they report to be underestimated by up to 
0.5-1 orders of magnitude, implying an even stronger disagreement with the galaxy-scale 
measurements. 
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Fig. 7. — SFR surface density Esfr vs. gas surface density S gas . The dashed line represents 
the Kennicutt- Sch midt relation, while the lower dotte d line represents the observational fit by 



Bigiel et al. 



( 120081 ) and the top dotted line is the fit by lWu et al. 



for indivi dual low- to intermediate -mass star-forming regions by 



bols) and 



(200 5tL We also plo t the data 



Evans et a 



Heiderman et al 



(120 lOl ) (x symbols) and for mass ive clumps by 



(120 lOl ) (* symbols). The filled square rep resents OMC-f (se e 



and the filled triangle is Taurus (see, e.g.. 



Heiderman et al. 



(2009 ) (+ sym- 



Heiderman 



Vazquez-Semadeni et al 



et al. 



20oJ 



20101 ) . The arrows in the massive 



clumps from the latter authors indicate the likely correction to the SFR due to their appli- 
cation of extragalactic methods to Galactic regions. The solid black line shows the evolution 
of our calibrated model with i? inf = 10 pc. 

These cloud-scale observations occupy a well-defined locus in E gaa — Esfr space, which 
can be compared with our model. For this task, we choose an inflow radius i?i n f = 10 pc, 



for which our mode 
median mass of the 



reaches a maximum mass of M max « 2000 M w , almost identical to the 



Evans et al. 



(120091 ) sample. In Fig. [7] we t hen p lot th e evolution of thi s 



model, as well as the loci of the clouds from the lEvans et al.l (120091 ) and 



Heiderman et al. 
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( 120 lOl ). adding an upwards-pointing arrow to the latter points, of length corresponding to 
one order of magnitude, to indicate the likely underestimation of t he SFR for massive clumps 



by th e latter authors. We al so plot the data 



Heiderman et al. 



rom O MC-1 (from 



Vazquez- Semadeni et al. 



20101 ) . The model evolves from low to high 



20091 ) and Taurus (see e.g. 
values of both S gas and Esfr- 

It is interesting to note that the model passes first through the locus of the low-mass 
star-forming clouds and later near the locus of the clumps forming massive stars. This 
means that the model predicts that present-day, relatively quiescent, low-mass-star forming 



clouds may evolve into massive-star-forming ones in a fe w to several Myr. A sim i 



conclusion was reached through numerical simulations by 



Vazquez-Semadeni et al. 



ar 



(120091 ) . 



This reinforces the idea that the dispersion of the obse rvational data is due to different 



evolutionary states of the clouds in a sample (see, e.g., 



Bigiel et al 



2010). 



5.4. Stellar age distribution 



One important prediction of our model is that the SFR increases over time. This is 
because, as the cloud contracts and its mean density increases, the fraction of star- forming 



gas in the c 



( 11999 



2000 



oud in creases. An increasing SFR has already been proposed by 



Palla fc Stahler 



20021 ) on the basis of the age distribution in various low- and high-mass 



clusters. However, this result has been questioned, since there is evidence suggesting that 



the older st ars are not genuine members of the clusters, but rather 



popu lation (jHartm ann 



20081 ). Moreover, 



2003 



3allesteros-Paredes fc Hartmann 



2007 



jelong to a different 



Heitsch fc Hartmann 



Hartmann! (120031 ) has posed the conundrum that, if most clouds form at 
an accelerated pace only over the last few Myr, and form stars at a very slow rate over the 
previous 10 Myr or so, most clouds should be found to be in the slow-star-forming period, 



but this is not what is observed. Our evolutionary scenario for clouds may offer a solution 
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to this debate. 

Because in our model we compute the total mass of stars (AM) formed at each time 
step, we can readily obtain the total number of stars formed during that time step as the 
integral over all masses of the IMF, normalized to AM. The left panel of Fig. [8] shows the 
stellar age histogram for our calibrated model, with i? inf = 10 pc (M max « 2OOOM ) at the 
end of its life - i.e., when it has completely lost its gas. We show the histograms for a case 
with feedback off (model LN-FO) and one with feedback on (model LN-F1). It is clear from 
this figure that indeed the age distribution is concentrated towards young ages, although a 
few older stars exist. 




Age [Myr] Age [Myr] 



Fig. 8. — Left: Stellar age distribution for our calibrated model with R ini = 10 pc (M max rs 
2000M q ), calculated at the end of the cl oud's evolution. Righ t: Compilation of the age 



histograms for the associations studied by 



Palla fc Stahleii f l2000l ) 



r 



These results can be compared with the age histograms presented by iPalla fc Stahler 



( 120001 ) for the Orion Nebula Cluster (ONC), the Taurus- Auriga region, Lupus, p-Oph, 
Chameleon, Upper Scorpius, and IC348. In these clusters, the fraction of stars with ages up 
to 1 Myr ranges from ~ 30% to 66%, while the fraction of stars with ages up to 4 Myr is in 
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the range 80 - 97%. Moreover, only in the case of Upper Scorpius does the age histogram 
peak at an age larger than 1 Myr, namely at 3 Myr. For this association, the fraction of 
stars with ages < 1 Myr is only 11%, while the fraction with ages < 4 Myr is 71%. We 
show a compilation of these in the right panel of Fig. |HJ 



We can see th at, qualitatively. 



;he st ellar age histogram at the end of our model's life 



resembles those of iPalla fc Stahlerl ( |2000| ). although, quantitatively, the model's histogram 
in Fig. [8] is much more concentrated towards short ages. However, this must be due to the 
fact that it was calculated at the end of the model's evo l ution . Clearly this is not the case 



for the clusters and groups analyzed by 



Palla &: Stahlerl ( 120001 ). because, as those authors 



themselves point out, in most cases the clusters are still embedded in their parent clouds, 
with only Upper Scorpius being already exposed. This means that we should consider our 
model before the end of its life. 

In Fig. [9] we show the age histogram for the calibrated model, calculated at 1 and 
2 Myr before the end of its evolution, and compare it with one of the histograms from 



Palla fc Stahlerl (120001 ) - that for p-Oph. We see that the histogram becomes less peaked 



as earlier times befo re the cloud's de s tructi on are taken, becoming more closely resemblant 
to the histograms of IPalla fc Stahlerl (120001 ) . This is because, as the SFR increases towards 
later times, the fraction of young stars becomes increasingly larger. In particular, the 
histograms for 1 and 2 Myr before the cloud's dispersal seem to bracket the histogram 
for p-Oph. As seen from the right panel of Fig. [HI the other regions are less concentrated 
towards short ages. According to our model, then, p Oph is a somewhat more evolved 
region, well matched by our model at < 2 Myr before its dispersal, while the other regions 
would correspond to somewhat less evolved stages, at slightly earlier times (3-4 Myr) before 
destruction. 



These results suggest a possible resolution of the debate between the Palla-Stahler and 
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Fig. 9. — Stellar age distribution for our calibrated model with i?i n f = 10 pc (M max rs 
2000M Q ), calculated at one and two Myr before the end of the cloud's evolution, com pared 



with the corresponding distribution for the p-Oph association (IPalla fc Stahler 



2000). 



the Hartmann groups. Specifically, although our model indeed predicts an inc rease of the 
SFR in collapsing clouds, this does not conflict with the conundrum posed by 



Hartmann 



(120031 ): no fully formed molecular clouds are observed without significant amounts of star 
formation because the clouds themselves are evolving. Thus, at the time when they had much 
lower SFRs, they were not fully formed yet, and thus not identifiable as large molecular 
clouds. Indeed, the clouds' mean density was lower, and thus, in reality, they probably 
consisted of a few molecular clumps immersed in a still-atomic interclump medium. Only 
in the last few Myr of their evolution, the clouds are dense enough on average that most of 
their bulk is already molecular, and by that time they are forming stars at a much higher 
rate, as observ ed. A similar conclusion has been recently reached on the basis of numerical 



simulations by 



Hartmann et al 
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6. Discussion 



6.1. The constant- .M rms assumption 



A feature of our model that may appear odd at first sight is that we have taken 
-Mmis = 3 as the (constant) fiducial value for the rms Mach n umber o f the t urbulence 



within the cloud, as it is contrary, for example, to the famous 



L arson! (Il98ll ) velocity 



dispersion-size scaling relation. However, it must be recalled that in this p aper we are 



specifi cally assuming that such a relation, or its more modern rendition by 



Heyer et al. 



( 120091 ). is a manife station of the gravitationa 



of the turbulence (IBallesteros-Paredes et al. 



contraction of the cloud, rather than a feature 



20111 ) . Thus, the relevant rms Mach number 



must be the remainder after the collapsing motions have been removed. A competition 



may be set up between the transfer of kinetic energy from t 



he co 



turbu lent ones and the dissipation ( iVazquez-Semadeni et al. 



1998 



l apsing motions to the 



Klessen Sz Hennebelle 



20101 ) . and so, in the absence of a reliable model, we consider t hat the assumption of a 



constant rms Mach number with the value typical for the CNM (jHeiles fc Troland 



20031) is 



reasonable, although a possible alt ernative recipe for it s initial value would be to take it 



equal to the inflow Mach number ( IBanerjee et al 



2009). We consider that further work is 



necessary to better constrain this parameter. 



6.2. The lognormal PDF assumption 

A similar situation arises for the density PDF in the cloud, which we have assumed to 



powe r -law tail at high densities 



2008 



Kainulainen et a 



However, 



Kritsuk et al. 



2009 



Klessen 



2000 



Kritsuk et al. 



Dib k Burkert 


2005; 


Vazauez-Semade) 


ai et al. 


2011 




Ballesteros-Paredes et al. 


2011b 


)• 



(120111 ) have suggested that such power-law tails are the effect of 
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the development of local collapsing sites with power-law density profiles. In this case, as 
explained in §2.5} the power-law tail in the PDF would be the result of the collapse, rather 
than the seed for it, and thus the relevant PDF for the seeds for future collapse should be 
the underlying lognormal one, after removal of the already-collapsing regions. Moreover, 
there is no complete theory for how the density PDF should evolve in time from a lognormal 
to a power-law. Unknowns such as the timescale for the transition, the density at which the 
power-law tail starts, and the final slope of this region are uncertain at present. As above, 
we consider that further work is necessary to clearly resolve this issue, and in the meantime 
we settle for the lognormal PDF assumption. 



6.3. Accelerating star formation 

An important precision is in order concerning the acceleration of star formation in 
our model. Indeed, our model predicts that the star formation rate increases over time, 
and therefore, star formation (SF; strictly speaking, the instantaneous stellar mass, Mg(t)) 
accelerates. However, it is common to find the statement in the literature that it is the 
SFR that accelerates. This is not the case for our model with feedback. The SFR is the 
time derivative of Mg. Since the SFR increases in time, the second time derivative of Mg is 
positive, and thus the SF accelerates. However, the second time derivative of the SFR (the 
third derivative of Mg) is negative for our model with feedback (see the left panel of Fig. H]), 
and thus strictly speaking the SFR decelerates. 



6.4. Room for improvement 



In the present model, we have bypassed the supporting effect of all forms of pressure, 
and replaced it by the empirical "Larson factor", /l, which effectively lengthens the 
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timescale for collapse. In the case of the original work by iLarsonl (119691 . Appendix C), this 
factor represented the support from thermal pressure which, incidentally, should be most 
im portant during the earlier stages of the collapse. Calibrating against the SAF1 simulation 



by 



Vazquez-Semadeni et al. 



(120101 ). we found a value of /l rouhly 8% larger than the 
one found by Larson, suggesting perhaps that turbulent pressure added a certain (small) 
amount of support (the magnetic field was not included in that simulation). Including 
physically-motivated terms into the collapse prescription that account for thermal, turbulent 
and magnetic support is an important goal, which we will attempt to pursue in a future 
contribution. 

Nevertheless, it is interesting that our model, calibrated in the non-magnetic case of the 
SAF1 simulation, gives a good match to a number of observational properties of molecular 
clouds in a wide range of masses. This suggests that magnetic support is not crucial in 
these objects. In turn, this is consistent with the recent realization that star-forming cloud s 



tend to be magnetically s upercritical in general (IBourke et al. 



Troland fc Crutcher 



2001 



Crutcher et al 



2003 



20081 ). and thus they should be essentially in a free-fall regime. 



Finally, in this paper we have not considered the effect of supernova explosions towards 
the late evolutionary stages of the clouds. This may help in reducing the model's SFR at 
those stages, probably bringing it to better agreement with the observations (cf. Fig. [7]). 



7. Summary and Conclusions 



by 



In this paper we have develop ed a semi-empirical analytical model (based on simulations 



Vazquez-Semadeni et al 



20101 ) in which a MC is formed by converging WNM flows. We 
assumed that the inflow collision produces a CNM cloud, through nonlinear triggering of 
the thermal instability, and that the cloud becomes turbulent through the combined action 
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of the latter and various other dynamical instabilities, such as the nonlinear thin shell, 
Kelvin-Helmholtz, and Rayleigh- Taylor ones. We assumed that the rms Mach number of 
this turbulence remains fixed at the typical values in the CNM (Ai s ~ 3), and that over 
its evolution, the cloud develops further nonthermal motions related to its collapse, n ot its 



internal turbulence. We also assumed that the cloud forms stars with a iKroupal ( l200ll )-type 
IMF, so that massive stars only appear when a sufficiently large number of stars has formed 
to adequately sample the high-mass tail of the IMF. Finally, we assumed that the density 
PDF in the cloud has a lognormal shape and a fixed width (corresponding to a constant 
turbulent Mach number M. Tms ), but whose maximum shifts towards higher densities as the 
cloud contracts and becomes denser on average. 

Using the same WNM infl ow parameters as the simulation labeled SAF1 from 



Vazquez-Semadeni et al. 



( 120101 ) . namely _R in f = 64 pc, n- mi = 1 cm , and M. x 



we calibrated the model by searching the density threshold for star formation, tt-sf, that 
best matched the simulation's evolution of the SFR and SFE. Our calibrated value was 
nsF = 10 6 cm -3 . With the nsF, Jiinf, and M. xms parameters fixed, the only remaining free 
parameter of the model is the WNM inflow radius i?i n f, which essentially controls the 
maximum mass reached by the model cloud, M max . Varying this parameter we then match 
the model to clouds of various masses, and compare with various properties of such clouds. 

The generic behavior of the model cloud, with the parameters of the SAF1 simulation, 
are as follows: i) The size of the model cloud decreases faster (by gravitational contraction) 
in a case without stellar feedback (model LN-FO) than in a case with it (model LN-F1). 
This is because feedback partially evaporates the cloud, thus reducing its gravitational 
potential, and slowing its collapse, ii) The model without feedback approaches a final state 
of zero size with finite mass (a singularity), while the case with feedback approaches a final 
state of zero mass at finite size; i.e., it is completely consumed by SF and ionization before 
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if reaches zero size, iii) Although the SFR increases in both cases, it accelerates over time 
in the LN-FO model, while it decelerates over time in the LN-F1 model. 

We then set out to apply the model to explain a number of observed features of 



molecular clouds of a wide range of masses. First, we compared the predictions o 



Kawamura et al. 



our model 



d2009|), in 



with the evolutionary scenario for GMCs recently proposed by 
which the GMCs start out having no massive stars, then have reduced numbers of them, 
so as to only have isolated HII regions, and finally have large numbers of them, so as to 
clearly contain massive clusters. We find that our model, with a value of i?i n f that gives 
M max ~ 10 5 M , comparable to the mass range reported by those authors, spends similar 
times in each of the stages reported by them. 

We also investigated a model cloud with R ini = 10 pc, corresponding to M max ~ 
2000 M . We find that such a model cloud evolves in the EsFR-S gas , or Kennicutt-Schmidt, 
diagram, in such a way that it passes first through the locus of individual low-mass-star 



for ming clouds and la ter t hrough the locus o 



by 



Evans et al 



(j2009h and 



Heiderman et al. 



high -mass-star forming clumps, as reported 



10) 



(120101 ) . Next, we compared an evolved stage of 



the calibrated model, also using R ini = 10 pc, with the physical conditions in the OMC-1 
massive clump, finding that it has similar physical conditions after ~ 26 Myr of evolution 
since its parent molecular cloud first formed, although it spends only about 2 Myr in a state 
comparable to OMC-1. 

Finally, we investigated the stellar age distribution in our isolated-cloud model with 
i?i n f = 10 pc, showing that, taken a few Myr before the en d of the cloud's l ife, it is 
consistent with the corresponding distributions presented by 
various clusters and associations. Furthermore, the model predicts that the shape of this 
age distribution depends on the evolutionary stage of the system, being more peaked 



Palla k Stahler! (120001 ) for 



towards young ages as the system grows older, because of its increasing SFR. 
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We conclude that our evolutionary and collapsing model of molecular clouds adequately 
represents actual clouds of a wide range of masses, with no need whatsoever for the 
consideration of equilibrium states. In this sense, the present model, although idealized, 
represents a promising first attempt at a non-equilibrium model for molecular clouds and 
their star-forming properties. 

We acknowledge helpful discussions with Javier Ballesteros-Paredes and Jose Franco. 
An anonymous referee suggested to perform the comparison against the SAF1 simulation, 
helping us to better constrain the model. This work has received partial finacial support 
from CONACYT grant 102488 
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